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Abstract 

We study the influence of a time-dependent potential on the motion of solitons in a quasi 
one-dimensional Bose-Einstein condensate by solving the corresponding Gross-Pitaevskii equa- 
tion. For a suitable choice of the external potentials as well as the initial soliton characteristics 
time-dependent trapping of the soliton in a prescribed subarea of the condensate can be achieved. 
Adiabatic perturbation theory is shown to work remarkably well for large switching on times of 
the trapping potential and allows to perform a detailed study of the degree of trapping in the 
complete phase space of the soliton center. A remarkable spiral pattern of the degree of trapping 
as a function of the soliton characteristics is observed and explained. 



PACS numbers: 03.75.Fi,05.45.Yv,05.30.Jp,02.30.Jr 



INTRODUCTION 



Since the experimental discovery of Bose-Einstein Condensation (BEC) in dilute alkali 
metal vapours in 1995 this field has experienced an enormous increase of interest (for 
recent reviews see refs.[0, |3|]). Indeed, the perspectives are very promising and range from 
coherent matter wave optics such as atom lasers, interferometry or atom lithography to 
precision measurements and quantum information processing. External static electric and 
magnetic or electromagnetic fields are of equal importance to the production as well as ma- 
nipulation of the condensed phase. Using different traps it is nowadays possible to produce 
in particular low i.e. (quasi-) one- and two-dimensional condensates M. 

One of the most promising ideas for the manipulation of coherent matter waves is the 
so-called atom chip which consists of a network of charge and current carrying elements, 
such as wires, on, for example, a semiconductor surface || (for a comprehensive review see 
R. Folman et al. 0). The electric and magnetic fields generated by the corresponding 
charge distributions and currents can then, together with homogeneous external fields, be 
employed to guide and control the coherent atomic ensemble which is moving a few micron 
above the surface. In particular the different field configurations offer the possibility to 
create microtraps, waveguides and other atom-optical devices on a single atomic chip. Very 
recently the first BEC on such an atom chip has been prepared J7|, |8], [5], |10[ using microscopic 
traps. Subsequently the condensate has been coupled to an atomic conveyer belt |nj which is 
used to transport the condensed cloud nondestructively over a macroscopic distance parallel 
to the chip surface J7J. 

According to the above optics with atoms being either in the thermal or condensate phase 
is now well-established. An important question posed here is whether one can do optics 
with excitations of the condensate. Specifically we will study the possibility to control the 
motion of nonlinear excitations of the condensate, in general, and in particular of solitons. 
Experimentally there exist several quantum-phase-engineering techniques to generate dark 
solitons in a Bose-Einstein condensate (see, e.g., ref. fl2fl ) which, according to theoretical 
predictions oscillates in the BEC. The question then arises how one could influence 
or even guide their motion, e.g. by modifying the above mentioned domains, so as to 
transfer "information" via solitons in a controllable manner. Here we provide a step in this 
direction by investigating the possibility of trapping the motion of a soliton in a prescribed 
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subarea of the condensate via a time-dependent potential. The time-dependent change of 
the potential can, for example, be implemented on atom chips by switching currents flowing 
in the corresponding wires. As we shall see the degree of trapping m t in the subarea, that is 
a measure for the strength of the confinement of the motion of the soliton, shows a twofold 
dependence. First it depends on parameters such as the switching time of the potential. 
Second it turns out that, for fixed parameters, the degree of trapping m t is a function of the 
phase space of the soliton center showing a beautiful spiral pattern. 



SOLITONS AND THE GROSS-PITAEVSKII EQUATION 

The above condensates consist of dilute ensembles of weakly interacting alkali atoms. 
In spite of its weakness the interaction is responsible for many of the properties of the 
condensate. In particular the diluteness at the nano-Kelvin scale implies that the dominant 
s-wave scattering of binary collisions provides a scattering length a which is much smaller 
than the mean interatomic spacing. As a consequence it is for many purposes a good 
approximation to introduce a contact potential as an effective interaction. In a mean- 
field description this yields the Gross-Pitaevskii equation (GPE) (see ref.0 and references 
therein) 

ih%-*(r,t)= (-^-V 2 + V ext + g\V(r,t)Ay(r,t) (1) 
ot \ 2m J 

where m is the mass of a single atom, g = is the coupling constant and V ext represents 
an external potential. In our case V ext will contain both the external time-independent 
potential which provides the overall confinement of the BEC as well as a time-dependent 
part to be specified below. 

Due to the nonlinearity occuring in eq. ([]]) the GPE is capable of describing energetically 
low-lying collective excitations of the condensate such as solitons. For repulsive atom-atom 
interaction (g > 0), which is the situation investigated here, only dark solitons occur. Soli- 
tons are well-known to play an important role in traditional nonlinear (fibre) optics. In 
the case of coherent matter waves the existence of solitons has been demonstrated recently 



followed by a number of investigations on their properties [|T3| , |14[ P~5f1 . 

Eq.([j]) provides a mean-field description of the condensate and its excitations in three 
space dimensions. However, present day experiments allow to prepare in particular low- 
dimensional condensates of (quasi-)two or one dimensional 0] character showing very inter- 



esting thermodynamical as well as microscopic properties. In the case when the confinement 
for two of the three spatial dimensions is much stronger than in the third dimension, the 
GP equation can be reduced to an effective quasi- (1 + l)-dimensional GP equation fll6| . 
For repulsive inter-atomic interactions, the latter equation can be expressed in the following 
dimensionless form, 

iu t + lu xx -\u\*u = VUx)«, (2) 
where the spatial coordinate x and time t are normalized to the harmonic oscillator length 



a± = yh/mu± and oscillation period, 1/lu± respectively. The frequency u± belongs to 
the two dimensions with strong confinement. The field u describes the macroscopic wave 
function \1/ of the condensate, according to the following scaling relation 

*<*■*> =(isir «<*■«> <3) 

The external potential V ext is decomposed into a time-independent part V\ which ensures 
the overall confinement for all times and an explicitly time-dependent part according to the 
following appearance 

V ext (x,t) = V^x) + f(t - t )V 2 (x) (4) 

The additional potential V 2 is turned on adiabatically. The major change of the potential 
occurs within the time interval [—t c /2 + 1 , t c /2 + t Q ] where t a is the characteristic switching 
on time and t c is the duration of the switching on process. Asymptotically (t — > oo) the 
full potential V\ + V 2 has been turned on. In the case under consideration the function f(t), 
which plays the role of a switch, is chosen to be 

f{t) = l( 1+Ua &(tzh)) (5) 



2 V ~ y t c 

In the limit t c — > f(t) leads to the step function 9{t) i.e. we encounter an instantaneous 



switching on process 



V 1 (x)+V 2 (x), t>t 



It is important to note that the role of the additional potential is, as we will see, to 
confine the motion of the solitons within a certain spatial region, which can be varied 
almost arbitrarily by the particular choice of V 2 (x). There are two possibilities concerning 
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the effect of V^(x) on the shape of the initial trapping potential Vi(x), namely to add either 
a barrier or a well to the confining potential. In the former case the spatial extension of 
the ground state wave function is not affected significantly; however, in the latter case, the 
spatial extension of the new ground state is clearly reduced leading, in the extreme cases of 
very deep and narrow wells, to the collapse of the condensate. We will therefore consider 
only the case where V%(x) adds a barrier to the confining potential V\(x). 



COMPUTATIONAL APPROACH 



The GPE in eq.([|) is solved numerically using the Split-Step Fourier method [17]]. In the 
numerical simulations we set up a solitary excitation for the initial time t = using the 
following Ansatz for the wave function of the soliton |L9[ 



u(x, 0) 



(cos (f)(0) tanh£ + i sin 0(0)) 



f = (x-x o (O))cos0(O) (7) 

where \x is the chemical potential of the BEC, x o (0) is the position of the soliton center 
at t = and 0(0) the corresponding soliton phase angle (|0(O)| < |). Assuming that the 
soliton trapping potential V 2 is switched on much later than t — 0, i.e. t Q ^> 0, we can take 
the function U(x, 0) in eq.(^) equal to the BEC trapping potential V\(x). The expression ([?J) 
is based on the results of adiabatic perturbation theory for dark solitons fl8 l, which can be 



applied in the case where V\(x) varies slowly on the soliton scale, and turns out to be a very 



good estimate for the initial soliton wave function [[19] . Solving eq. (g) numerically, we then 
obtain the evolution of the soliton on top of the condensate. The fact that the initial form 
of the soliton (0) is an approximation to the exact soliton wave function introduces, besides 
the standard numerical propagation error, an additional source of inaccuracy with respect 
to the soliton dynamics. However, as can be seen in the following, the time evolution of the 
Ansatz (^) is very close to the expected solitary wave propagation. If the switching on time 
interval t c is large then the potential varies slowly with time and the adiabatic perturbation 
theory can be used to estimate, at least approximately, the soliton + BEC dynamics. In fact 
we will see later that the critical value for t c below which the adiabatic perturbation theory 
breaks down is t c pa 10. Particularly, the equations describing, within adiabatic perturbation 
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theory (for t c > 10), the motion of the soliton center [JT3, [HJ are given by: 

d 2 x Q dU 2^JI 



in 



m eff = —?T7^ ( 8 ) 



dt 2 dx D ' 11 COS 2 0(0) 

and are asymptotically valid, i.e. for t — > — oo or t — > oo where the time dependence of the 
external potential can be ignored. For slowly varying switching on function f(t) (t c > 10) 
one can take for the effective potential U = V ex t where V ex t is the total external potential. 
Equations (§) can, in general, be solved numerically to obtain the trajectory of the center 
of the soliton. The results obtained by the direct numerical integration of the GPE are 
in a fairly good agreement with the ones obtained by the numerical integration of eqs.(||). 
Therefore, the latter can safely serve as a guide for the choice of the appropriate initial 
conditions for potential trapping used in the much more expensive numerical propagation 
of the GPE. 

Having fixed the form of the external potential the only free parameters entering in the 
solution of eq.fl2|) are the initial values x (t = 0), (pit = 0) of the position of the center and 
the phase of the soliton, respectively. After switching on the potential V^x), the possible 
trapping of the dark soliton in a prescribed region of the condensate depends strongly on 
a; o (O),0(O). The amplitude, and therefore also the velocity, of the soliton are determined 
through sin (f)(0). Large values for the amplitude lead to a slowly moving soliton. If the 
soliton is moving very slowly the time needed to perform a full oscillation within the con- 
densate is comparable with the lifetime of the condensate itself: In this case trapping is not 
an issue. If the amplitude of the soliton is small, the velocity is large and the trapping of 
the soliton becomes impossible. One could think of avoiding the loss of trapping by using 
a more pronounced external potential. However this leads in general to a highly unstable 
condensate wave function. As a consequence there are several competing factors which de- 
termine the evolution of the composite system (soliton + condensate) and its properties: (a) 
the initial conditions of the solitonic excitation (x o (0), 0(0)) controlling the evolution of the 
soliton, (b) the form of the potential determining the shape of the condensate and (c) the 
duration t c of the switching on process of the soliton trapping potential V^- If the potential 
is turned on suddenly the BEC becomes very unstable and there is no initial configuration 
for the soliton leading to trapping before the BEC is destroyed. Therefore we expect to 
have trapping of the soliton, which can be clearly observed within the BEC lifetime, only for 
restricted values of x a , 0(0) and t c as well as special forms of the soliton trapping potential 
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Vo. 



GENERAL ASPECTS OF TRAPPING 

Before proceeding with the discussion of our numerical results let us define the measure 
m t for the degree of trapping for a given trajectory of the center of a soliton in a subarea of 
the BEC 

AH - AM 

mt = (9) 

where A^~\ AH are the amplitudes of the oscillation of the center of the soliton for times 
t <C £ and t t , respectively. We will investigate soliton trapping for two different cases of 
the external potential V ext . The condensate trapping potentials Vi(x) possesses in both cases 
a single minimum. After turning on the soliton trapping potential the total potential V ext 
will exhibit two minima for the external potential V^} and correspondingly three minima 
for V^l ( see fig s - 3 and 5, respectively). In order to reduce the number of free parameters in 
our study we fix the values for the parameters of the external potential such that the BEC 
ground state matches well the dimensions of the experimentally prepared condensates. 
Let us first consider the external potential 

f) = ax 4 + f(t - t )([3x 2 + 7 ) (10) 

where 7 = 0.2, f3 = -0.001, a = 1.25 • 10" 6 , t c = 40 and t = 600. Having fixed the 
external potential we remain with the free parameters x o (0) and 0(0) which determine the 
solitonic excitation. First we calculate the degree of trapping m t for a 100 x 100 grid in the 
(x o (0), 0(0)) space. The results of these calculations are shown in the contour plot of Fig. 1. 
They are obtained employing adiabatic perturbation theory. For a subset of the gridpoints 
we have also calculated, as a test, m t using the full GPE propagation. The results turn out 
to differ at most by 20% for this value of t c (The integration of the GPE is computationally 
too expensive in order to obtain a well-resolved sin(<p ) — x Q plot in Fig.l). We observe 
a spiral pattern which determines the region with significant trapping (m t > 40%) of the 
soliton in the condensate. There are regions with trapping of the order of 100% depicted 
by the white areas in the arms of the spiral pattern. We have repeated our calculations 
for larger as well as for smaller t c . For increasing t c (we have chosen t c = 70) the arms 
of the spiral become thicker and finally they meet each other forming a connected almost 



rectangular region where trapping occurs. Decreasing t c (t c = 15) we observe that the arms 
of the spiral shrink more and more and the trapping region tends to disappear. 

The origin of this spiral pattern of the trapping measure m t can be understood qualita- 
tively by recalling that the switching on process is accompanied by an energy transfer to the 
soliton. The soliton center trajectory moves for t <C t , i.e. before the trapping potential is 
turned on (V — Vi), on an closed phase curve in the phase plane (x a , sin 2 (p). During the 
switching on process this phase curve is deformed. In case of strong trapping for t ^> t 
(V — Vi + V 2 ) the resulting phase curve is a closed phase curve again but with a significantly 
smaller size. In other words the available phase space volume for the motion of the soliton 
center has shrinken significantly due to the trapping. The transition between these two 
phase curves occurs via an inward spiral which represents the natural transient between the 
two closed curves of different size. Quantitatively the spiral pattern can be understood by 
calculating the energy transfer in e.g. the limit of a sudden switching on and imposing the 
final condition of a strong confinement i.e. trapping. 



REFLECTION AND TRAPPING OF SOLITONS 



The basic mechanism leading to the soliton trapping is the reflection of the soliton at 
the walls of the external potential. At the reflection point of the soliton center its kinetic 
energy vanishes and the soliton center therefore possesses only potential energy. During 
this reflection process the phase of the soliton is changed dramatically. Such a reflection 
process is monitored in Fig. 2a where three different time instants (before, at and after the 
reflection on the central hump of the external potential) are depicted in the trajectory of 
the soliton center obtained by numerical integration of the GPE. In Fig. 2(b,c,d) we show 
the phase function of the soliton for these three time instants, as well as the corresponding 
probability density \u\ 2 . In particular, for t = 982 (before the reflection), Fig. 2b shows the 
soliton phase as a function of x, which is characterized by a smooth jump from a region of 
low phases (left of the soliton center) to a region of high phases (right of the soliton center). 
As the soliton approaches the barrier its depth increases and the phase function becomes 
steeper. As a result, at the reflection point shown in Fig. 2c (t = 1041), the phase jump 
becomes step-like, i.e. discontinuous and after the reflection (see Fig. 2d for t = 1104) the 
regions of low phases and high phases are interchanged while the change of the phase across 
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the soliton center becomes smooth again. These three snapshots are calculated using the 
same parameters for the external potential as in Fig. 1. 

In Fig. 3a we show snapshots of the soliton+BEC dynamics for twelve subsequent time 
instants using the initial conditions (x o (0), sin 0(0)) = (—5, V0.15) and t c = 40, t a = 600. 
The soliton trapping is clearly visible for times t > 600 for which the motion exclusively takes 
place in the right well of the potential. The small fluctuations due to radiation of the BEC 
can be well distinguished from the solitonic excitation up to times t w 1200. For significantly 
larger times the size of the fluctuations becomes of the order of the soliton amplitude finally 
leading to the decay of the condensate. A rough estimation of the trapping degree in this 
case gives m t « 0.6. This can be also seen in Fig. 3b where the corresponding trajectory 
of the soliton center (solid line) is shown. In the same plot we present also the results of 
adiabatic perturbation theory for a comparison (dashed line). It is readily seen that the 
results of perturbation theory fit very well the numerical calculations. The relative error 
in the present case in the amplitude of the soliton oscillations is less than 30%. We have 
repeated our calculations for two other choices of the switching on parameter t c . The results 
of the calculations are given in Fig. 4(a,b). In Fig. 4a we show the trajectory of the soliton 
center for the potential ( |T0"D using t c = 15. Presented are both the GPE integration results 
(solid line) as well as the adiabatic perturbation theory results (dashed line). It is clearly 
seen that in this case (small t c ) the discrepancy between perturbation theory and exact 
numerical integration is much larger than for t c = 40 in figure 3. The opposite behaviour is 
observed in Fig. 4b where we present analogous results for the case t c = 70. 

The above shape of the trapping potential does not allow for the possibility of trapping 
at the center of the BEC. Taking into account that laboratory BECs possess a finite lifetime, 
it is easier to experimentally preserve a soliton at the center of the BEC where boundary 
effects are minimal. Therefore we have studied the effect of soliton trapping for the potential 

V(x, t) = ax 6 + f(t- t ) ((3x 4 + jx 2 ) (11) 

To obtain experimentally accessible length and time scales for the BECs size and lifetime 
we have used the parameter values 7 = 0.00398, (3 = -1.27510" 5 , a = 1.0210~ 8 , t c = 100 
and t Q = 950. After turning on the potential V2 the resulting total potential possesses three 
minima. However the BEC lifetime turns out to be very sensitive with respect to both the 
turning on time t c and the initial properties of the solitonic excitation. To achieve trapping of 
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the solitonic excitation in the central well of the total external potential one has to start with 
a black soliton (i.e. a soliton with zero initial velocity). Even small deviations of the initial 
velocity of the soliton of the order of 10~ 2 are sufficient to prevent trapping. In Fig. 5(a,b) 
we show the evolution of a solitonic excitation with initial values (x o (0), sin 0(0)) = (—10, 0) 
and the external potential given by eq. (JTT|) . Fig. 5a shows two snapshots describing the 
inital state as well as the final state of the composite system. We observe the trapping of 
the initial soliton in the central well of the potential for times t > 1000. BEC radiation 
effects are visible but can be considered as small perturbations for times t < 1500. For 
larger times these perturbations increase substantially finally leading to the destruction of 
the condensate. The corresponding trajectory of the soliton center (solid line) is given in 
Fig. 5b together with the results of adiabatic perturbation theory. A very good agreement 
is also here observed. This is expected due to the large value of t c (t c = 100). The trapping 
degree is m t « 0.56. 

Some additional comments regarding the stability of the BEC are in order. The numerical 
simulations reveal that the stability properties of the BEC depend strongly on the value for 
t c in the switching on function /(£). For small values of t c the system becomes very unstable 
while large values of t c a significant stability is achieved. For t c ~ 10 we encounter a 
practically stable BEC+soliton system. 

CONCLUSIONS 

In summary, we have demonstrated that a switching on process of a suitably chosen 
external potential in a quasi one-dimensional BEC leads to trapping of solitonic excitations 
in a prescribed subarea of the condensate. Depending on the initial soliton's position and 
velocity as well as the shape of the external potential the trapping i.e. the confinement of 
the motion of the soliton, can be very strong (degree of trapping m t ~ 1.0). For an external 
potential with two minima there is a finite region of initial conditions in the phase space 
of the soliton center for which trapping occurs. If the switching on process is slow enough 
(i.e. of the same order as the period of the soliton oscillations) then the corresponding 
composite system (BEC+excitations) is stable for long time intervals and time-independent 
perturbation thoery predicts satisfactory the soliton dynamics. The degree of trapping shows 
a spiral pattern in the phase space of the soliton which can be understood in terms of the 
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energy transfer during the switching on process. Using an external potential with three 
minima the trapping of the soliton becomes a much more difficult task. Only solitons which 
are initially close to black ones can be trapped in this case; this can be understood by the 
fact that shallower solitons have enough initial kinetic energy to jump the barrier induced by 
Vi- The above results suggest that it is possible to monitor and control solitonic excitations 
in BECs by means of proper matter wave devices. 
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FIGURES 



FIG. 1: The contour plot of the trapping degree mt as a function of the phase space coordinates 
(x o (0), sin 2 0(0)) of the initial solitonic excitation for the external potential (10) and t c = 40. The 
calculations are made within adiabatic perturbation theory. 



FIG. 2: Illuminating the soliton reflection process: (a) Three time instants before (t = 982), at 
(t = 1041) and after (t = 1104) the reflection of the soliton at the barrier of the external potential 
[see eq.([lO|)], depicted at the trajectory of the soliton center. The parameters used for the numerical 
integration of the GPE are the same as in Fig. 1, for a dark soliton with initial position x Q = — 5 
and velocity sin 0(0) = v0T5. (b) The density |u(x,i)| 2 and the phase function arg(u) for t = 982 
(before the reflection), (c) The same as in (b), but for t = 1041 (at the reflection). Here the density 
of the soliton becomes minimum and the phase jump across the soliton becomes step-like, (d) The 
same as in (b), but for t = 1104 (after the reflection). The regions of low and high phases have 
been interchanged. 



FIG. 3: (a) The density |u(x,t)| 2 for twelve different times for the potential in eq. © and t c = 
40. The initial conditions for the solitonic excitation are (x Q (O),sm0(O)) = (-5,v / 0T5). (b) 
The corresponding trajectory of the center of the solitonic excitation. Both the numerical GPE 
integration results (solid line) as well as the adiabatic perturbation theory results (dashed line) are 
presented. 
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FIG. 4: (a) The trajectory of the center of the solitonic excitation for the external potential ( |l~0| ) and 
t c = 15. The solid line corresponds to the numerical integration of the GPE while the dashed line 
is obtained using adiabatic perturbation theory. The initial conditions for the solitonic excitation 
are (x o (0), sin(/)(0)) = (—5, V0.15). (b) The same as Fig. 4a with t c = 70. 



FIG. 5: (a) The density |ti(x,i)| 2 for two different times t = 200 and t = 1000 for the potential 



in eq.(ll). The initial conditions for the solitonic excitation are (x o (0), sin(f>(0)) = (—10,0). (b) 
The corresponding trajectory of the center of the solitonic excitation. Here also the solid line 
corresponds to the results obtained by numerical integration of the GPE while the dashed line are 
the adiabatic perturbation theory results. 
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